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1.0 SUMMARY 


Program L216 (DUBFLX) is structured as nine overlays, one main and eight primary 
overlays. Input into the program is made via card and magnetic files (tapes or disks). 
Output from the program consists of printed results and magnetic files containing 
geometric and aerodynamic data. 

Although L216 serves as a module of a DYLOFLEX system, it can be operated as a 
standalone program. Subroutines used by L216 include routines embedded in the 
program code and routines obtained from the standard FORTRAN and DYLOFLEX 
libraries. 



2.0 INTRODUCTION 


Program L216 (DUBFLX) can be used either as a standalone program or as a module of 
a program system called DYLOFLEX (see fig. 1), which was developed for NASA under 
contract NASl-13918 (ref. 1). DUBFLX is used to compute subsonic unsteady 
aerodynamic forces which are required in aeroelastic analyses as developed under the 
criteria given in reference 2. The current version is a derivative of the original Air 
Force H7WC program (ref. 3). Structure of program H7WC was modified and 
streamlined for NASA. Program L216 contains enhancements to the streamlined 

version, as well as revisions to meet the DYLOFLEX specifications^ in order to make it 
compatible with the other modules of the DYLOFLEX system. 

This volume contains information for maintaining and/or revising the program L216. To 
meet this objective, the following items are described in some detail. 

• Program design and structure 

• Overlay purpose and description 

• Input, output, and internal data base descriptions 

• Test cases 


^Clemmons, R. E.: Programming Specifications for Modules of the Dynamic Loads 
System to Interface With FLEXSTAB. NASA contract NASl-13918, BCS-G0701, 
September 1975. (Internal document) 
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Figure 1. — DYLOFLEX Flow Chart 











3.0 PROGRAM DESIGN AND STRUCTURE 


Program L216 (DUBFLX) has the overlay structure shown in figure 2. The main 
overlay, (L216,0,0), calls into execution the primary overlays required to perform the 
analysis defined by card input data. To solve the non-AIC problem, the following 
primary overlays are called: 


(L216,l,0) 

INITIAL 

(L216,2,0) 

PRT2 

(L216,4,0) 

GEQ 

(L216,5,0) 

AMSOL 

(L216,7,0) 

GIN 

(L216,8,0) 

SAMODE 

To solve the AIC problem, th< 

(L216,l,0) 

INITIAL 

(L216,2,0) 

PRT2 

(L216,3,0) 

GEWSSL 

(L216,5,0) 

AMSOL 

(L216,6,0) 

AICSUB 


Each of the overlays is described separately in the following sections. Included are the 
overlay’s purpose, file communication, macro flow chart, and a list of subroutines 
called. 

The external files read/written by L216 are displayed in figure 3. 

Program L216 can serve as a module of the DYLOFLEX system or as a standalone 
program. When the program is run by itself, it becomes the user’s responsibility to 
generate input data in the format required by L216. (See vol. I of this document for the 
required data and formats.) 

This program requires subroutines that are not part of the L216 code. Some routines 
are automatically obtained from the standard FORTRAN library when the program is 
loaded. Others, however, are stored in the DYLOFLEX alternate subroutine library, 
which must be declared at the time of loading. Subsequent sections describe each 
overlay separately and contain tables displaying the routines called and the library in 
which they are located. All subroutines have only one entry point. 

This volume describes the program in a macro sense. A more detailed discussion 
appears in the comments contained in the program source code. Each routine contains 
a preface describing . the routine’s purpose, author, analytical steps, modification 
history, input data, output data, glossary of variables, and list of other routines called. 
Embedded within the executable code are comments labelling each section and 
explaining logical branches. 
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( From L215 (INTERP). 
the Modal Interpolation 
Program (ref.^ ) 


Figure 3. — External Data Files Communicating With L216 (DUBFLX) 


6 






3.1 OVERLAY (L216,0,0) - L216vc 


The main overlay of L216 is named L216vc. The last two characters of the program 
name define the program version and the last correction set against that version, 
where v is a letter indicating the program version and c is an integer number indicating 
the last correction against version v. 

Purpose of L216vc 

L216vc performs certain bookkeeping tasks, directs the execution of the primary 
overlays, and aids the communication between primary overlays via labelled common 
blocks. 


Analytical Steps of L216vc 

L216vc performs its task in the following steps: 

1. L216vc calls the subroutine PRGBEG to place the program header on the printed 
output. 

2. The sequential and random access files used as scratch files are defined to the 
program. 

3. The overlay INITIAL reads card input data for geometry and modal data. 

4. For modal data saved on the file SATAP, the overlay SAMODE interpolates for 
mode shapes at all aerodynamic elements. 

L210VC repeats steps 5 through 8 for each reduced frequency: 

5. L210VC computes the downwash factors matrix [D] in the overlay PRT2. 

6. It computes the modified kinematic downwash, the [W^ic] matrix for AIC (in 
overlay GEWSSL), or the non-AIC [w] (in overlay GEQ) problems. 

7. It solves the [D^] [W] = [ACp] equation in overlay AMSOL to obtain the 
pressure coefficients. 

8. It prints aerodynamic results for AIC problems by calling overlay AICSUB. For 
non-AIC problems, it calls overlay GIN to print the pressure and generalized air 
forces and sectional derivatives. 

9. It looks for a CASE card to initiate the analysis of a new data set. 

10. It calls the subroutine PRGEND to place the program termination header on the 
printed output. 

The macro flow chart of this overlay appears in figure 4. The subroutines called are 

displayed in table 1. 
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Figure 4. —IMacro Flow Chart of Overlay (L216,0,0)—L216vc 
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Figure 4.— (Concluded) 
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Table 1.— Routines Called by OVERLAY (L216,0,0,0)—L216vc 


OVEPLAY (L216, 0,0) 

PROGRAM L216VC 

ELK DAT 
CLOSMS* 

FBTAD+ 

FETDEL+ 

FILLM 

IRDCPD+ 

OPENMS* 

PRGBEG* 

PRGENC+ 

RETURNF* 

SECOND* 

L216VC ’ KRTETP+ 

INITIAL OVERLAY (L21 6, 1,0) 

PRT2 OVERLAY (L21 6, 2,0) 

GEWSSL OVERLAY (L21 6, 3,0) 

GEQ OVERLAY (L216, 4,0) 

AMSOL OVERLAY (L21 6, 5,0) 

AICSUB OVERLAY (L2 16, 6,0) 

GIN OVERLAY (L21 6, 7,0) 

SAMODE OVERLAY (L 2 16, 8,0) 

+ indicates a routine in the DYLOFLEX library 
* indicates a routine in the FORTRAN library 
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I/O Devices of L216vc 


L216vc reads from the card input file card sets numbered 1.0 through 18.0. L216vc 
writes on the printed output file, messages indicating the program name and version, 
as it begins and ends execution. All other input/output to and from L216 takes place in 
lower level overlays. 


3.2 OVERLAY(L216,1,0) - PROGRAM INITIAL 


Purpose of INITIAL 

INITIAL reads the card input data defining the problem geometry and modal data. It 
stores geometry in labelled common blocks for later use in L216 and writes geometry 
matrices onto the file NGETP for subsequent use in the DYLOFLEX equations of 
motion (L217 EOM) program described in reference 5, 


Blocks of card input data are introduced by the keywords: 


REDUCED FREQUENCIES To introduce k-values 


GEOMETRY 

PANEL 


To introduce geometry data 

To define the lifting surface panel geometry and modal 
data 


BODY 

MODE 


To define the slender body geometry and modal data 
To specify the modal data set 


Figure 5 is a macro flow chart of the program INITIAL. Table 2 displays the 
subroutines and functions called by the program INITIAL. 


I/O Devices of INITIAL 


INITIAL reads from; 

INPUT 
(card reader) 

INITIAL writes on: 

OUTPUT 

(printer) 

NGETP 

(user-named file) 

NTP8 and NTP9 
(random access 
scratch files) 


Card sets 4 through 17.4 


The interpreted card input data and error 
diagnostics 

The geometry data in matrix form 


The labelled common blocks of geometry and modal 
data are saved temporarily on these files to lower 
core storage requirements 
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Figure 5. —Macro Flow Chart of Overlay (L216, 1 ,0}— INITIAL 
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Table 2.— Routines Called by OVERLAY (L216,1 ,0)— INITIAL 


I 

\ 


\ 

\ 

i 

i 

i 


?! 
6 ' 


OVERLAY (L2 16^1,0) 
PROGRAM INITIAL 


FETAE+ 

FETDEL+ 

{ WRITMS* 
STINDX* 
GESTORE WRTETP+ 
NAMFIL+ 


POLYB 


INITIAL < 


PRT1 


PEDB 


STINDX* 
< WRITMS* 
READMS* 
FILLM 
' READMS* 
- WRITMS* 
STINDX* 

V 


WRITMS* 


* indicates a routine in the DYLOFLEX litrary 

♦ indicates a routine in the FORTRAN library 
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3.3 OVERLAY(L216,2,0) - PRT2 


Purpose of PRT2 

PRT2 computes the steady and oscillatory components of the matrix of downwash 
factors [D] for all the lifting surface elements. It also modifies the downwash factors for 
the two planes of symmetry (Y = 0 and Z = 0). The matrix is written on scratch tape 
NTPl for later use. 


Analytical Steps of PRT2 

PRT2 computes the downwash factors matrix from the expression 



y-7}, Z-f, Cl), M) do 


sending 

element 


where K is a nonplanar kernel function and cr is the tangential spanwise coordinate. 
The term ai is the frequency of oscillation and M is Mach number. 

Figure 6 is the macro flow chart of PRT2. Table 3 displays the subroutines called by 
PRT2. 


I/O Devices of PRT2 


PRT2 reads from: 
NTP8 

(random access 
scratch file) 

PRT2 writes on: 

OUTPUT 

(printer) 

NTPl 
(sequential 
scratch file) 


The labelled common block data 


Error diagnostics and intermediate debugging data 
if requested 

The [D] matrix 


3.4 OVERLAY(L216,3,0) - GEWSSL 
Purpose of GEWSSL 

GEWSSL computes the downwash column vectors {WaIc} for each of the preselected 
AIC modes. Figure 7 is the macro flow chart of GEWSSL. 

GEWSSL calls the subroutine READMS (in the FORTRAN library) and SSL. 
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Figure 6. —Macro Flow Chart of Overlay ( L21 6,2,0)— PRT2 
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1 


PFT2 

♦ indicates 

* indicates 


Table 3— Routines Called by OVER LA Y ( L21 6,2,0}— PRT2 


INCRO 


KERNEL 
* IDF1 
IDF2 


" SNPDF 
REALMS* 

REQFL* 

a routine in the DYLOFLEX library 
a routine in the FORTRAN library 


16 




Figure 7. —Macro Flow Chart of Overlay ( L21 6,3,0)— GEWSSL 
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I/O Devices of GEWSSL 


GEWSSL reads from: 


NTP8 

(random access 
scratch file) 

GEWSSL writes on: 

OUTPUT 

(printer) 

NTP3 
(sequential 
scratch file) 

NTP4 
(sequential 
scratch file) 


The geometry data GEWSSL writes on: 


The modal data, the augmented downwash factors 
matrix, and error diagnostics 

The downwash matrix [W^jc] for all strips 

The downwash matrix [Waic] columns 
3.5 OV&RLAY(L216,4,0) - GEQ 


Purpose of GEQ 

GEQ computes the downwash matrix [W] for all non-AIC analyses. For the structure 
with slender bodies, it also computes the pressure coefficients [ACp] and the induced 
downwash [Wb] on the lifting surfaces. The modified downwash matrix [W-Wb] is 
then computed and stored. 


Figure 8 is a macro flow chart of GEQ. Table 4 displays the routines called by GEQ. 

I/O Devices of GEQ 


GEQ reads from: 


NTP3 The modified downwash matrix [W-Wb] 

(sequential 

scratch file) 


NTP8 

(random access 
scratch file) 

GEQ writes on: 
OUTPUT 

NTP4 
(sequential 
scratch file) 
NAETP 

(user named file) 


Labelled common block data 


The calculated results and error diagnostics 

The modified downwash matrix [W-Wb] by 
columns 

The induced downwash matrix [Wb] 
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Figure 8. —Macro Flow Chart of Overlay ( L21 6,4,0)— GEQ 












Table 4.— Routines Called by OVER LA Y (L216,4,0)—GEQ 


GEQ 4 


♦ indicates 

* indicates 


CWGBD <1 


BPBNT 
COMB IN 
READMS* 

SUMER i TKER 

I WRTETP+ 

STINDX 


AUGP 


f READMS* 
4 SUMER 


BPRNT 


► 

< 


PWGBD i 


TKER 

WRTETP+ 


I COMBIN 
y READMS* 

READMS* 

routine in the DYLOFLEX library 
routine in the FORTRAN library 
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3.6 OVERLAY(L216,5,0) - AMSOL 


Purpose of AMSOL 


AMSOL solves the simultaneous equations represented by the augmented 
factors matrix, [D/W] It forms the [D'^W] matrix, which is used to obtain 
pressure coefficients [ACp] . 


downwash 
the lifting 


Analytical Steps of AMSOL 

AMSOL checks the core storage available, initializes scratch files, and calls the 
subroutine QUAS for the matrix solution 

where [D] is the downwash factors matrix and [W] is the modified downwash matrix. 
Figure 9 is the macro flow chart of AMSOL. Table 5 displays the subroutines called by 
AMSOL. 


[ACp] = [D]-l [W] 


I/O Devices of AMSOL 


AMSOL reads from; 


NTPl 
(sequential 
scratch file) 

The matrix [D] 

NTP4 
(sequential 
scratch file 

The matrix [W] 

AMSOL writes on: 


OUTPUT 

(printer) 

The AIC solution (if requested) 
diagnostics 

NTP3 
(sequential 
scratch file) 

The matrix [D"' W] 

NAETP 
(sequential 
output file) 

The quasi inverse matrix [D]’^ 


AMSOL writes/reads partitions of the [D/W] matrix on/from these sequential scratch 
files: 
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Figure 9. —Macro Flow Chart of Overlay ( L21 6,5,0)— AMSOL 








Table S.-Routines Called by OVER LA Y (L216,5,0)-AMSOL 


AMSOL 


FETAD<- 

FETDEL+ 

- QUAS 

REQFL+ 
RETURN F+ 


indicates a 

routine 

in 

the 

DYLOFLEX 

library 

indicates a 

routine 

in 

the 

FORTRAN 

library 



NTP2 

NTPIO 

NTPll 

LTAPE 

3.7 OVERLAY(L216,6,0) - AICSUB 
Purpose of AICSUB 

AICSUB calculates and prints the AIC stability derivatives and pressure coefficients. 
When requested, it also computes and prints the harmonic gust coefficients. 

Figure 10 is the macro flow chart of AICSUB. Table 6 displays the routines called by 
AICSUB. 

I/O Devices of AICSUB. 

AICSUB reads from: 

NTP3 The solution matrix [D'^W] for analysis 

(sequential 

scratch file) 

NTP8 The geometry data and modal data (integration 

matrix B) 

NTP9 The saved AIC dimensional matrix to form the AIC 

force coefficients 

AICSUB writes on: 

OUTPUT The results and error diagnostics 

(printer) 

NTP9 The AIC dimensional matrix (this data is 

subsequently read by subroutine FIN) 

3.8 OVERLAY(L216,7,0) - GIN 

Purpose of GIN 

GIN computes and prints the pressure coefficients, sectional and total lift and moment 
coefficients, and generalized air forces. It also writes on the file NAETP the 
aerodynamic results required by the DYLOFLEX equations of motion (L217 EOM) 
program (ref. 6). 

Figure 11 is the macro flow chart of GIN. Table 7 displays the subroutines called 
by GIN. 
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Figure 10. —Macro Flow Chart of Overlay (L216,6,0)—AICSUB 
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Table S.-Subroutines Called by OVERLA Y ( L21 6,6,0}— AlCSUB 


AICSOE 


♦ indicates 


AEB 


AI 


STINDX* 
' REAOMS* 
WPITMS* 

L 

STINDX* 
- HEADMS* 
WR ITMS* 


FIN 

PEADMS* 


routine in the FORTRAN library 
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Figure 11. — Macro Flow Chart of Overlay (L216,7,0j—GIN 
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Table 7.— Routines Called by OVERLA Y (L216,7,0)—GIN 

f 

REAEWS* 

STINDX* 

GIN ^ WRITMS* 

WRTETP* 

■k 

+ indicates a routine in the DYLOFLEX library 
♦ indicates a routine in the FORTRAN library 
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I/O Devices of GIN 


GIN reads from: 
NTP3 
NTP7 
NTP8 
NTP9 

GIN writes on; 

OUTPUT 

(printer) 

NAETP 
(sequential 
output file) 


The solution matrix [D ^W] 

The slender body pressure coefficients 
The geometry data 
The modal data 

The calculated results and error diagnostics 

The calculated aerodynamic data required by the 
equations of motion (L217 EOM) program (ref. 6) 

3.9 OVERLAY(L216,8,0) - SAMODE 


Purpose of SAMODE 

SAMODE computes the modal input data from the SA arrays generated by the 
DYLOFLEX interpolation program (L215 INTERP) (ref. 4 ). For each lifting surface or 
slender body SAMODE reads the SA array from the file NTPSA and calculates the 
displacements, slopes, and the integration matrix [BQ] . 

Figure 12 is the macro flow chart of SAMODE. Table 8 displays the subroutines called 
by SAMODE. 

I/O Devices of SAMODE 


SAMODE reads from: 
NTPSA 


The modal interpolation array 
{SA} from L215 (INTERP) (ref. 5) 


NTPl 
(sequential 
scratch file) 

SAMODE writes on: 

OUTPUT 

(printer) 


The interpolated modal data 


The interpolated modal data and 
error diagnostics 
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ENTRY 
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Figure 12. —Macro Flow Chart of Overby (L216,8,0)—SAMODE 



















Table 8.— Routines Called by OVER LA Y ( L21 6,8,0)— SAMODE 


SAMODE i 


EXPMD 

FETAD+ 

FETDEL+ 


MESA 


AINTG+ 

< DEE I VS 
READTP+ 


SEADMS* 

PEADTP+ 

EETURNF*- 

STINDX* 

WPITMS* 


indicates a 

routine 

in 

the 

DYLOFLEX 

library 

indicates a 

routine 

in 

the 

FOETPAN 

library 



NTP8 The modal displacements and 

(random slopes 

access file) 


The integration matrix [BQ] 


3.10 DATA BASES 

L216 data bases include input and output files plus internal scratch files and common 
block storage. 

3.10.1 INPUT DATA 

The input data are from two sources: cards and magnetic files. The geometry data is 
always read from cards. The modal data may be input on cards or on the magnetic file 
generated by the DYLOFLEX modal interpolation (L215 INTERP) program (ref. 4). 

Card Input Data 

The card input for L216 is described in detail in section 6.5.1 in volume 1 of this 
document. 


NTP9 
(random 
access file) 


Magnetic File Input Data 

L216 reads only one file generated externally, the file NTPSA, containing modal 
interpolation data generated by the program L215 (INTERP). L215 and the file are 
described in detail in reference ^ . 

3.10.2 OUTPUT DATA 

The results of an L216 execution will be printed and, optionally, written onto magnetic 
files. 


Printed Output Data 

A sample of the printed output of L216 is included with the sample problem in section 

7.2 of volume I. 


Magnetic File Output Data 

Two magnetic files, NGETP and NAETP, are generated by the program. NGETP 
contains the saved geometry data, and NAETP contains the aerodynamic data. A 
detailed description of formats and contents of both appears in section 6.4 of volume I of 
this document. 
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3.10.3 INTERNAL DATA 


Two methods are used to transfer data from one portion of the program to another; 
common blocks and scratch (temporary) storage files. In some parts of the program, 
even the common block data are stored temporarily on the scratch files to reduce the 
memory required to run the program. 

Common Blocks 

The common blocks used inside L216 may be grouped into four types: 

1. Some labelled common blocks are defined in the (0,0) overlay and always reside in 
core. These labelled common blocks are named: 

/APPLE / 

/CHKPRT / 

/CNTRL / 

/INOUT / 

/MAE / 

/NTPS / 

/PROB / 

/PRTFLG / 

/PRTNQ / 

/RAF / 

2. Several other labelled common blocks are defined in more than one overlay, but not 
in the (0,0) overlay. These common blocks are stored temporarily on the random 
access scratch file named NTP8 and read into core for the overlays in which they 
are required. They are named: 

/BODY/ 

/GEBW/ 

/PTIQ/ 

/xxz/ 

/XYZ/ 

/YZY/ 
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3. 


A few labelled common blocks are defined in only one overlay, and are used to 
communicate with and between subroutines of that overlay. The locally defined 
common blocks are: 


/MGE / 

/PIGW / 

/RWBUFF/ 

/STUFF / 

4. Blank common is defined in several primary overlays. It is used to store large 
arrays. The array sizes are variable, dependent upon the problem sizes. 

Table 9 shows which common blocks are used in each overlay. Table 10 contains 
detailed descriptions of the labelled common block elements. The following 
abbreviations are used in table 10 to describe the elements of each common block. The 
column labelled "T” in table 10 indicates variable type. The possible codes and their 
meaning in a FORTRAN sense are: 

R Real 
I Integer 
C Complex 
L Logical 
H Hollerith 

If the dimension column is left blank, "1” should be assumed. 
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Table 9. —Common Blocks Defined in Each Overlay 


COMMON 

BLOCKS 


OVERLAY ^ 


(L216,0,0)-L216vc 


(L216,2,0)-PRT2 


(L216,3,0)-GEWSSL 

(L216,4,0)-GEQ 

(L216,5,0)-AMSOL 


(L216,6,0)-AICSUB 

(L216,7,0)-GIN 

(L216,8,0)-SAMODE 


/APPLE/ 

/BODY/ 

/CNTRL/ 

/CHKPRT/ 

/GEBW/ 

/MAE/ 

/MGE/ 

/NTPS/ 

/INOUT/ 

/PROS/ 

X X 

X XXX 

X X X X X 

X X X X X 

X XXX 

XXX 

XXX 

XXX 

X X X X X 

XXX 


XXX 

X X 

XXX 

X X X X X 

XXX 

X XX 

X X 


O \ 

a \ \ 

2: o IS \ 
Eh Eh O P4 

« « Eh T- < 

ft K 

\ \ \ \ \ 


X 

X 

XXX 

X 

X 

X 


X X X X 


X X 


/STUFF/ 

/XXZ/ 

/XYZ/ 

/YXY/ 





















Table 10.— Contents of Common Blocks 


LABELED COMMON NAME: APPLE 


DESCRIPTION: 

Constants commonly used. 















ENG. 


NO. 

VARIABLE 

T 

DIMENSION 

NOM. 

DESCRIPTION 

1 



PI 

R 


TT 

Constant 

2 

TWOPI 

R 


2tt 

Constant 

3 

PI02 

R 


tt/2 

Constant 

4 

PISQ 

R 


2 

TT 

Constant 

5 

RAD EG 

R 



Conversion factor (radians to 






degrees ) 






,, -4 

6 

SMALL 

R 


€ 

Constant (10 ) 


36 











Table 10.— (Continued) 








Table 10.— (Continued) 










Table 10.— (Continued) 


LABELED COMMON NAME: CNTRL 


DESCRIPTION: Problem control flags and counters 


NO, 

VARIABLE 

T 

DIMENSION 

ENG. 

NOM. 

DESCRIPTION 

1 

NDELT 

I 


6 

Symmetry flag 

2 

IQ 

I 



Loop counter 

3 

IR 

I 



Loop counter 

4 

JSPECS 

I 



Symmetry Flag (Plane Z=0) 

5 

NGUST 

I 



Gust Flag 

6 

NYAW 

I 



Stability Derivative override 
flag . 

7 

NAIC 

I 



AIC Analysis Flag 

8 

NMD IN 

I 



Type of mode input 

9 

NPFM 

I 



Flag to print aero results 

10 

JRF 

I 



Reduced frequency index 

’ 3 

IMX 

I 



Number of AIC modes 
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Table 10.— (Continued) 


LABELED COMMON NAME: 


GEBW 



NO. VARIABLE 


DIMENSION NOM. 


DESCRIPTION 
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Table 10.— (Continued) 


LABELED COMMON NAME: INPUT 


DESCRIPTION: Identifiers and Error Counter 


NO. 

VARIABLE 

T 

DIMENSION 

ENG. 

NOM. 

DESCRIPTION 

1 

INFIL 

I 



INPUT file (card reader) 
TAPES 

2 

lUTFIL 

I 



OUTPUT file (line printer) 
TAPE6 

3 

lAERO 

I 



Flag to 
results 

save aerodynamic 
on the file NAETP. 






=0 do not save results 






=1 save generalized forces and 
unsteady pressures 






=2 save 
[D)“ 

the matrices 
^ and [FI 






=3 save 
both 

all items listed under 
1 and 2 above. 

4 

KERROR 

I 



Fatal error counter 

5 

ICASE 

I 



Current 

case number 

6 

ICOND 

I 



Current 

condition number. 
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Table 10.— (Continued) 


LABELED COMMON NAME: MAE 


DESCRIPTION: 

Control matrix for saved 

Aerodynamic results 










NO. 

VARIABLE 

T 

DIMENSION 

ENG. 

NOM. 

1 

DESCRIPTION 

■ 

MAC 

I 

85 


See Section 3.6 of Volume I 
for details. 

2 

SHFTC 

R 

3,40 


Shift for panel coordinates for 
modal interpolation. 









Table 10.— (Continued) 


LABELED COMMON NAME; MGE 


DESCRIPTION: Geometry data to be saved on NGETP 


NO. I VARIABLE I T | DIMENSION] NOM. 


DESCRIPTION 


1 MGC 


2 MPSG 


3 MISG 


4 MSBG 


Control matrix (see Section 3.6 
of Volume I for details) . 

Geometric data for primary lift- 
ing surfaces. 

Geometric data for interference 
lifting surfaces. 

Geometric data for slender 
bodies . 


5 NLIFS 


Number of primary lifting 
surfaces. 


6 NINTS 


Number of interference lifting 
surfaces . 


7 NSBDS I 


Number of slender bodies. 


8 MSBC 


I 400 


Box-surface correspondence 
array . 
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Table 10.— (Continued) 


LABELED COMMON NAME: NTPS 


DESCRIPTION: File names 


NO. VARIABLE 


DIMENSION 


DESCRIPTION 




8 NTPS 


File name SCTAPl 


File name SCTAP12 


File name SCTAP3 


File name SCTAP4 


File name SCTAP5 
File name SCTAP6 
File name SCTAPl 3 
File name SCTAP8 


File name SCTAP9 


10 NTPIO 


11 NAETP 


12 NGETP 


File name SCTAPIO 

User-supplied file name of 
saved aerodynamic results. 

User-supplied file name of 
saved geometry data. 


13 NTPQ 


14 NTPSA 


Not used. 

User-supplied file naime of the 
modal SA arrays. 
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Table 10.— (Continued) 


I 


LABELED COMMON NAME: PROB 

DESCRIPTION: 

General Data 

Array 















ENG. 


NO. 

VARIABLE 

T 

DIMENSION 

NOM. 

DESCRIPTION 

1 

NBOX 

I 



Number of lifting surface boxes. 

2 

NB 

I 



Number of slender bodies. 

3 

NBE 

I 



Number of slender body elements. 

4 

NTOT 

I 



Total number of aerodynamic 
elements (NBOX+NBE) . 

5 

NMD 

I 



Number of modes. 

6 

NRF 

I 



Number of reduced frequencies. 

7 

NSTRIP 

I 



Number of strips on lifting 
surfaces. 

8 

NSV 

I 



Number of strips on plane Y=0. 

9 

NBV 

I 



Number of boxes on vertical 
surfaces in plane Y=0. 

10 

NBF 

I 



First box number on vertical 
surfaces in plane Y=0. 

11 

NBL 

I 



Last box number on the vertical 
surfaces in plane Y=0. 

12 

NPC 

I 



Alternative mode selector of 
AIC analysis. 

13 

NMTP 

I 



Total number of modal polynomial 
coefficients defined on lifting 
surface panels. 

14 

NMTB 

I 



Total number of modal polynomial 
coefficients defined on slender 






bodies . 

15 

FMACH 

R 


M=v/c 

Mach number. 
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Table 10.— (Continued) 


LABELED COMMON NAME; PHOB 

(Conf d) 


DESCRIPTION; 

















NO. 

VARIABLE 

T 

DIMENSION 

ENG. 

NOM. 

DESCRIPTION 

16 

BETA 

R 


g = yi_M2 

Mach number function. 

17 

ACAP 

R 


A ref 

Reference Area. 

18 

BR 

R 


b 

Reference length used in comput- 
ing reduced frequency. 

19 

AKDBR 

R 


k/b 

"Normalized" reduced frequency. 

20 

RFREQ 

R 

70 

k=ub/V 

Reduced Frequency. 

21 

REFCHD 

R 


c 

Reference chord length. 

22 

REFSPN 

R 


s 

Reference semispan 

23 

NCARAY 

I 

40 


Niomber of chord-wise boundaries 
on a panel. 

24 

NSARAY 

I 

40 


Number of span-wise boundaries 
on a panel. 

25 

NBARAY 

I 

40 


Number of boxes on a panel. 

26 

WG 

R 



Gust velocity at XZRO. 

27 

VEL 

R 


V 

03 

Aircraft speed. 

28 

GZRO 

R 


rg 

Dihedral angle of gust refer- 
ence plane. 

2 9 

XZRO 

R 


XO 

Gust reference x location. 

30 

LIM 

I 

50, 3 


Summation limits of chord-wise 
integration on strips. 

31 

NSU 

I 



Number of lifting surfaces. 

32 

NSUBL 

I 

5, 80 


Lifting surface data. 
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Table 10.— (Continued) 


LABELED COMMON NAME: PRTFLG 


DESCRIPTION: Printing Options and Flags 


NO. 

VARIABLE 

T 

DIMENSION 

ENG. 

NOM. 

DESCRIPTION 

1 

NDATA 

I 



Flag to print downwash factors 
matrix [D] . 

2 

NSB 

I 



Flag to print slender body 
element downwash factor matrix 
[AZY] . 

3 

NPRl 

I 



Flag to print AIC analysis 
matrix solution [D“lw] . 

4 

NPR2 

I 



Flag to print AIC dimensional 
matrix [A] . 

5 

NPR3 

I 



Flag to print stability deriva- 
tives for preselected modes of 
AIC analysis. 
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Table 10.— (Continued) 



DESCRIPTION: Panel Definition Array 
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Table 10.— (Continued) 




ENG. 
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Table 10.— (Continued) 
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Table 10.— (Continued) 


LABELED COMMON NAME: RAF 


DESCRIPTION: Indices for Random Access Scratch Files 


NO. 


VARIABLE 


T 


DIMENSION 


ENG. 

NOM. 


1 

2 

3 

4 


IRAF8 

IRAF9 

FETGE 

BUFGE 


I 

I 

R 

R 


71 

71 

34 

514 


DESCRIPTION 

Index for file NTP8. 

Index for file NTP9. 

FET area for Saved Geometry file. 

Buffer area for Saved Geometry 
file. 
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Table 10.— (Continued) 


LABELED COMMON NAME; STUFF 

DESCRIPTION; Interprogram data used in finding downwash factors 


NO. VARIABLE 


DIMENSION 


DESCRIPTION 




n - 5 


y - n 


z - c, 


Y-PV See conunon blocks 

/GEBW/ and /PTlQ/ for 
Z-ZV . definition of these 
variables. 


Z-ZZl 



Z-ZZ2 


11 GAMSIG 


Dihedral angle of the receiving 
point. 

Dihedral angle of the sending 
point. 


Sin(Yr) Sine function of GAMS. 


Cos (Yr ) 


Cosine function of GAMS. 


17 10 


18 NBXS 


Sin(Yg) Sine of GAMSIG. 
Cos(Yg) Cosine of GAMSIG. 


Left-hand side/right-hand side 
of wing flag. 

Index 1 to number of boxes on 
the strip. 

Sum of sending point box numbers 
on wing. 
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Table 10.— (Continued) 


LABELED COMMON NAME: STUFF 

(Cont ' d) 



DESCRIPTION: 





















ENG. 



NO. 

VARIABLE 

T 

DIMENSION 

NOM. 


DESCRIPTION 

19 

NCPNB 

I 



Sum of sending point boxes on 
the panel. 

20 

NDBLE 

I 



Symmetry flag for plane Z=0. 

21 

SL 

R 



Sine of sweep angle. 

22 

CL 

R 



Cosine of sweep angle. 

23 

TL 

R 



Tangent of sweep angle. 

24 

TWOE 

R 


2e 

Width of sending box. 

25 

EE 

R 


e 

Semi-width of sending box. 

26 

E2 

R 


e2 

Square of width of sending box. 

27 

XMULT 

R 



Multiplication factor 
(SDELX/BPI) 

28 

ETA 

R 


y 

iKernel interpolation variables 

29 

ZET 

R 


Z 

) 


30 

ETASQ 

R 


-2 

y 

Square of ETA. 

31 

ZETSQ 

R 


-2 

z 

Square of ZET. 

32 

PARN 

R 





33 

ARE 

R 





34 

AIM 

R 




Coefficients of the parabola 

35 

BRE 

R 




. for planar part of kernel 







integration. 

36 

BIM 

R 





37 

CRE 

R 





38 

CIM 

R 
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Table 10.— (Continued) 







Table 10.— (Continued) 


LABELED COMMON NAME: XYZ 


DESCRIPTION: Strip Geometry Data 


NO. 


VARIABLE 


DIMENSION 


ENG. 

NOM. 


DESCRIPTION 


YS 

DELYS 

ZS 

DELZS 

FGAMMA 


R 

R 

R 

R 

R 


70 

70 

70 

70 

70 


6 


CWIG 


R 


70 


Y location of strip centerline. 
Strip location in Y plane. 

Z location of strip centerline. 
Strip width in Z plane. 

Strip dihedral angle. 

Average chord length of strip. 
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Table 10.— (Concluded) 
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Magnetic Files 


L216 uses two random access files (NTP8 and NTP9) and eight sequential files (NTPl, 
NTP2, NTP3, NTP4, NTP7, NTPIO, NTPll, and LTAPE) for temporary storage of data. 
The CDC random access routines READMS and WRITMS are used for NTP8 and NTP9. 
The sequential files are read/written with standard FORTRAN binary READ/WRITE 
statements. 


The files NTP2, NTPIO, NTPll, and LTAPE are used for the out-of-core solution of the 
problem data in OVERLAY(L216,5,0) - AMSOL. The files are continually written, 
rewound, read, rewound, rewritten, etc., during the solution of the [D/W] matrix. 

The file structure of the other scratch files is briefly described on the following pages. 

NTPl Scratch File 


When the modal data are defined by the SA array, SAMODE temporarily stores the 
interpolated modal data arrays on NTPl. They are reread in the same program, 
SAMODE. 

Usually, this file contains the downwash factors matrix [D] by rows. Tbe matrix is 
formed in the program PRT2. 


Matrix size 


(2*NBOX X l) 
(2*NBOX X l) 


Matrix 


{AWW}i 

{AWW}2 


( Real and imaginary 
values of [D] for the 
first row (first box) 


(2*NBOX X l) 


{awwInbox 


where NBOX = the number of 
box elements 



NTP3 Scratch File 


The contents of this file are dependent upon the user options selected. 

For AIC problems, NTP3 initially contains the gust modes and the downwash matrix 
[Waic]> created in GEWSSL and read by the subroutine SSL. 


Matrix size 
400 X 1 

300 X 1 
300 X 1 


300 X 1 


Matrix 



Gust righthand sides 

AIC downwash matrix 
(Mode 1) 

(Mode 2) 


(Mode NMD) 

where NMD = no. of modes 


After the AIC solution, NTP3 contains the AIC column of gammas formed by subroutine 
HSHDR. 


Matrix size 

(2*NBOx)x 1 


(2*NB0X)X 1 


Matrix 



Column 1 of gammas for AIC 
Column 2 


where NBOX - 
NMD = 


no. of box 
elements 

no. of modes 
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For non-AIC problems, NTP3 initially contains the modified downwash coliimn vectors 
created in subroutines AUGP or DWGBD and copied to NTP4 in the subroutine 
COMBIN. 


Matrix size 


(2*NB0X)X 1 
(2*NBOx)x 1 


(2*NBOX)x 1 


Matrix 



Modified downwash matrix 
(Mode 1) 

(Mode 2) 


(Mode NMD) 

where NBOX = no. of box 
elements 

NMD = no. of modes 


After the matrix solution, subroutine QUAS stores the pressure coefficients by columns 
for each mode. 


Matrix size 


(2*NBOx)x 1 
(2*NBOx)x 1 


(2*NB0X)X 1 


Matrix 



Pressure coefficients 
(Mode 1) 

(Mode 2) 


(Mode NMD) 

NBOX = no. of box elements 
NMD = no. of modes 
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NTP4 Scratch FUe 


NPT4 contains the modified downwash matrix [W-Wg]. The matrix is copied from 
NTP3 with no format changes by the subroutine COMBIN. 


NTP7 Scratch File 

This file contains the slender body pressure coefficients stored by the subroutine 
BPRNT. This file is read in the program GIN. 


Matrix size 


Matrix 


(2*NBE X)l 


(2*NBE)x 1 



Column of pressure 
coefBcients for all 
slender body elements 
for modes 1 through NMD 


where NBE = number of slender 
body elements 

NMD = number of modes 


File Structure of NTP8 

NTP8 is a random access file accessed through two index tables, IRAF8 and INDX8. The 
data written on this file is shown in table 11. The upper part of the figure displays the 
records defining polynomial modal data. The lower part displays records defining modal 
data input via SA arrays or tables. Table 12 displays the routines in which the data is 
read and written. 


File Structure of NTP9 

NTP9 is a random access file using two index tables, IRAF9 and INDX9. The contents of 
file NTP9 is shown in table 13. Table 14 displays the routines in which the data is read 
and written. 
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Table 11.— Contents of the File NTP8 


Array 

Size 

Indices 
in table 
IRAF8 

Indices 
in table 
INDX8 

ARP AY CONTENTS 

1000 

1 


Polynomial Modal data for Panels 

1000 

2 


Polynomial Modal data for Bodies 

950 

3 


AIC Modal Data 

420 

11 


Common Block /XYZ/ 

2 40 

12 


Common Block /BODY/ 

470 

13 


Common Block /XXZ/ 

3270 

14 


Common Blcck /GEBfc/ 

2400 

15 


Common Block /PT1C/ 

2 50 

16 


Common Block /YZY/ 

140 

17 


Common Block /PPTNQ/ 

1200 


1 

Displacements and Slopes for Mode 1 

1200 


2 

• 

Displacements and Slopes for Mode 2 

• 

1200 


• 

• 

NMD 

• 

• 

Displaceirent s and Slopes for Mode NMD 

NMD4-1 

1 


INDX8 table for Tabulated Mode 




Solution 
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Table 1 1.— (Concluded) 




Array 

Size 

Indices 
in table 
IRAF8 

Indices 
in table 
INDX8 

ARRAY CONTENTS 

1600 


1 

AIC Super columns for Strip 1 

1600 


2 

• 

AIC Super columns for Strip 2 

• 

1600 


• 

NSTRIP 

• 

AIC Super columns for strip NSTRIP 

NSTFIP+1 

5 


Random Index Table INDX8 

144 

4 


Static Stability Derivatives 
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Table 12.— (Concluded ^ 


Routine 

Performinq 

Read/Write 



Index 

Table 

Record 

Number 

R=Read 

w=write 



■ 



REP 

IRAF8 

■ 

3 

R 



■ 

4 

W 



■ 

5 

R 




s 



INDX8 

■ 

1 

R 



■ 

2 

R 




• 

• 



1 

• 

• 

• 

• 



■ 

. NSTRIP 

R 

OVERLAY (L216, 1 , 0) 










GIN 

IRAF8 


11 

R 




12 

R 




13 

R 




14 

R 




V 



where NMD 


= no. cf irodes 


N STRIP 


no 


of strips 














Table 13.— Contents of the File NTP9 


■ 1 

Array 

Indices 

Indices 

Array Contents 

Size 

in table 

in table 



IRAF9 

INDX9 





AIC Solution 

NSTFIP+ 1 

1 


INDX9 Table 

300 


1 

Integration Matrix Strip 1 

300 


2 

Integration Matrix Strip 2 

• 

• 

• 

300 


• 

• 

• 

NSTFIP 

• 

• 

• 

Integration Matrix NSTRTP 




Polynomial or Tabular Solution 

NMD + 1 

1 


INDX9 Table 

NTOT 


1 

Integration Matrix Model 1 

NTOT 

• 


2 

• 

Integration Matrix Model 2 

• 

• 

NTOT 


• 

• 

NMD 

• 

• 

Integration Matrix Model NMD 


NSTFIP 

= 

number 

of 

panel 

strips 

NMD 

= 

number 

of 

modes 


NTOT 

= 

number 

of 

boxes 

and body elements 
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Table 14.— Read/Write Activity on the File NTP9 



Routine 

Index 

Record 

R=Read 

Performinq 

Table 

Number 

W=Write 

Read/Write 




OVERLAY (L216, 1,0) 




GEE 

INDX9 

1 

W 



2 

w 



• 

• 

• 

• 



• 

• 



NSTRIF 

w 


IRAF9 

1 

w 

POLYE 

INDX9 

1 

R 



2 

R 



• 

• 



• 

• 



NMD 

R 


IRAF9 

1 

N 

RECE 

INDX9 

1 

W 

or 


2 

W 

OVERLAY (L216, 8,0) 


• 

• 

• 

• 

• 

• 

SAMODE 






NMD 

w 


IRAF9 

1 

u 
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Table 14.— (Concluded) 


Routine 

Index 

Record 

R=Read 

Performing 

Table 

Number 

W=Write 

Read/Write 




OVERLAY (L216,6,0) 

IRAF9 

1 

R 

AI 

INDX9 

1 

R 



2 

• 

R 

• 



• 

• 

NSTRIP 

• 

• 

R 

OVERLAY (L216, 7,0) 




GIN 

IFAF9 

1 

R 


INDX9 

1 

R 



2 

• 

R 

• 



• 

• 

NMD 

• 

• 

P 
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4.0 EXTENT OF CHECKOUT 


Brief descriptions of the test cases and the program options tested are listed below. The 

modal data for all test cases, except case 7, consists of two rigid body modes, pitching 

and plunging. 

Case 1: Original program H7WC sample problem consists of an aircraft fuselage as 

slender body, wing, and tail as lifting surfaces. Motion occurs in the Z 
direction only and is symmetric. 

Case 2: Two wing panels with Y symmetry and Z antisymmetry. Two reduced 

frequencies of 0.0 and 3.0 and Mach number of 0.130 are specified. 

Case 3: Two wing panels, a closed body with Z doublets. Two reduced frequencies of 

0.0 and 0.150 and Mach number of 0.130 are specified. 

Case 4: Fuselage, a nacelle, two wing panels, and a strut. Nacelle has Y and Z 

motions. Mach number is 0.85 and the reduced frequency is 0.50. 

Case 5: Two panel wings with AIC’s computed. Gust is also defined. 


Case 6: Two panel wings with tabular modal data. Same as test case 2. 

Case 7: Fuselage slender body with wing, horizontal tail, and vertical tail, as lifting 

surfaces. Motion occurs in the Y direction only, and is antisymmetric. The 
three modes are side translation, pitch, and yaw. 


Boeing Commercial Airplane Company 
P.O. Box 3707 

Seattle, Washington 98124 
May 1977 
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